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Abstract 

The effects of long- and stiort-range correlations on the charge distributions of some medium and 
heavy nuclei are investigated. The long-range correlations are treated within the Random Phase 
Approximation framework and the short-range correlations with a model inspired to the Correlation 
Basis Function theory. The two type of correlations produce effects of the same order of magnitude. 
A comparison with the empirical charge distribution difference between ^"^Pb and ^^^Tl shows the 
need of including both correlations to obtain a good description of the data. 



PACS numbers: 21.10.Ft, 21.60.-n 



1 Introduction 



The high precision of modern electron scattering experiments imposes severe constraints to nuclear models 
and theories. The study of the charge density distributions of medium and heavy nuclei by means of 
elastic electron scattering has allowed for a detail test of the capacities of the independent particle model 
(IPM) to describe the nuclear ground state Q. 

In the IPM the motion of each nucleon inside the nucleus is not affected by the presence of the other 
nucleons. The corrections beyond this picture take the generic name of correlations. The sources of 
correlations can be generically classified in two categories. The so called short-range correlations (SRC) 
are generated by the hard core repulsion of the nucleon-nucleon interaction which prohibits two nucleons 
to come too close to each other. The long-range correlations (LRC) are produced by collective excitation 
modes of the nucleus originated by the part of the nuclear hamiltonian neglected by the IPM. 

This separation in long- and short-range correlation is an artifact of our approach to nuclear structure 
usually starting from an IPM basis. Exact solutions of the many-body Schrodinger equation, ^ take 
automatically into account both kinds of correlations independently from their classification. When ap- 
proximate or effective theories are used, the role of one or the other kind of correlation can be emphasized. 
The results recently obtained with Fermi Hypernetted Chain (FHNC) techniques in doubly magic nuclei 
shows the lack of some correlation effects ^, |^, probably LRC. The aim of this article is to estimate 
the size of the effects of both kind of correlations on the charge density distributions of some medium 
and heavy nuclei. 

We estimate the LRC within a Random Phase Approximation (RPA) approach, and the SRC effects 
by using a first order expansion model based on the Correlated Basis Function (CBF) theory. We analyze 
both the quantitative and qualitative differences between the effects of the two types of correlations. The 
inclusion of both of them is necessary to describe the charge distribution difference between ^°^Pb and 
205rpj The main message of our work is that, in medium heavy nuclei, the two kind of correlations 
produce equally important effects beyond the IPM. 

The plan of the work is the following one. We briefly present in sects. || and |^ the theoretical 
background of our treatment of the correlations. In sect. ^ we first make a detailed description of the 
inputs of our calculations, then we present our results on the charge distributions of ^^O, ""^Ca and ^"^Pb. 
At the end of the section we show our results regarding the charge distribution differences between the 
above mentioned nuclei and their isotone partners with one proton less. Finally we calculate the charge 
density difference between the open shell nuclei ^"^Pb and 205^j g^j^j compare it with the empirical 
one. In sect, q we summarize our work and draw our conclusions. 



2 Long Range Correlations 

The RPA many-body states are defined as: 

QNl'foiRPA) 
and the operator Qn satisfies the equations of motions 
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for all the possible variations 6Qn- In the above equations we have indicated with H the nuclear 
hamiltonian, with [, ] the commutator and with lun the excitation energy. The RPA theory uses: 
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where and a are the creation and destruction operators respectively. The Xph{N) and Ypii{N) ampli- 
tudes are real numbers and are obtained by solving the RPA equations. 
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We are interested in evaluating the mean value of one-body operators on the RPA ground state defined 
by eq. In a second quantization language this can be written as: 



< O >=< -^oiRPAjlY^ < n\0\n' > alaf,,\-^o{RPA) > 



(5) 



where we have indicated with > a single particle (s.p.) wave function generated by the mean field 
basis. The matrix element < 'ifo{RPA)\aj^a^'\'ifo{RPA) > is evaluated as described in references 0, ^ 
and we obtain the result: 



< O < h\0\h > 



1 -JEE ivwi 



N p 



J2<p\o\p> 



^EEivwi^ 



N 



(6) 



where we have indicated with p and h the particle and hole states respectively. The operator defining 
the occupation number of a s.p. state fi is: 



^fu = (5a., 



(7) 



where a = p, h in the sums of eq. (g) . The modifications of the occupation numbers produced by RPA 
correlations has been studied in reference In this reference also the contribution of higher order 
excitations has been evaluated by means of a Second RPA calculation. 
We are interested in the charge operator: 



Q(r) = ^5(r-r,)i[l + T3« 



(8) 



with T3{i) = ±1 if the particle i is a proton or a neutron respectively. 

All our calculations are done in a j-j coupled spherical basis. This means that the s.p. wave functions 
are expressed as: 

,(r) = Rii,{r)Y^ < l^^ls\Jm. > Yi.^{?)Xs (9) 



where we have indicated with < | > the Clebsch-Gordan coefficients, with Yi^fj_(r) the spherical harmonics, 
with Xs the spin wave function and with t the third component of the isospin. 

Using the angular momentum coupled RPA equations and the expression for the s.p. states, we 
obtain the following expression for the charge density distribution of a doubly-closed shell nucleus: 
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where p and h indicate particle and hole states respectively limiting the range of the sums. In our 
calculations the sums on p and h are restricted only to proton states. We neglect the contribution of the 
neutrons to the charge distribution. In the IPM the Yph{J, N) amplitudes are zero and eq. ( p^ gives the 
well known expression of pipM- 

In this article we calculate also the charge distributions of nuclei with one proton less than the doubly 
magic ones. This is done within Landau-Migdal theory which allows to calculate the difference between 
expectation values of operators in doubly magic nuclei and nuclei with one nucleon less (or more) . The 
formalism has been widely used to calculate magnetic moments of nuclei around ^ospj-, g^j^j j^- j^g^g 
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Figure 1: Cluster diagrams contributing to the one-body density. The dashed hues represent the 
dynamical correlations, /i^, and the oriented lines the statistical correlations, po(ri,rj). A black dot 
associated with a point implies integration over its coordinates. 

been extended in references ^ to the calculation of the charge distributions. The basic equation to 
be used to evaluate this difference is: 



In the above equation V indicates the residual interaction used in RPA calculations, e the s.p. energy, 
ljJn the RPA excitation energy of the N-th state and \i > and \k > the s.p. states characterizing the 
transition in the A-1 system within the IPM. To calculate the difference between ground state charge 
distribution we consider J=0 and the fact that \i > = \k > is the last occupied proton s.p. level just 
below the Fermi surface. 

3 Short Range Correlations 

We evaluate the effects of SRC on the ground state charge distribution by using the approach of references 
1^ , The starting point is the following ansatz on the ground state: 



where we have indicated with $0 the Slater determinant of s.p. states. The correlation operator is written 
as a symmetrized product of two-body operators: 
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The expression of Fij used in the most sophisticated calculations, both in nuclear matter and in finite 
nuclear systems, Fij has the form: 

= E ^(^^3)0!, (14) 

p=l,6 

where 

Of=^'^ = [1, rr, ■ cr,, ] ® [1, • r,] (15) 

In the above equations we have indicated with Sij — (Sfy • cr^ f^j • cTj — cr.j ■ crj) the tensor operator. 
The expectation value of a generic operator O is: 

< O >= lli^li^ (16) 

In CBF theory the ground state is found through the minimization of the hamiltonian expectation value. 
The search for this minimum involves variations on the parameters of the correlation functions and of 
the mean field potential defining the s.p. basis. 

The charge distribution is obtained by inserting in eq. ( p^ ) the charge density operator (^) . We found 
convenient to rewrite the correlation function ( p^ ) as: 



G(l,2...yl) - 1 + 5[| 
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where we have defined — ~ 1 and h^-^^ = f^-^^. The expectation value is evaluated with standard 
cluster expansions techniques which allow us to express it as infinite sum of linked diagrams. The 
unlinked terms of the numerators are canceled by the denominator. 

At this point we introduce an approximation. We evaluate only those diagrams containing a single 
correlation function h^. The diagrams considered are shown in figure |l|. It is evident from Eq. (|l^) that 
the first term of our calculation is the IPM density which is already properly normalized. This implies 
that in the normalization procedure the contribution of all the correlation diagrams should vanish. It 
is shown in reference that the set of diagrams of figure |l] provides the correct normalization of the 
density. 

The validity of our approximation in calculating the charge distribution has been studied in reference 
Ip^ by comparing the results of the model with those of a Fermi Hypernetted Chain calculation. A test 
of the approximation in the evaluation of the nuclear matter charge responses has been done in reference 
Ip^ . More recently |^ the approximation has been tested in the evaluation of the one-body densities and 
the natural orbits. The conclusion of these investigations is that the approximation is reliable. This is 
probably related to the fact that the charge operator is a relatively simple one-body operator, therefore 
not very sensitive to the fine details of the many-body wave function. 

The calculation of the correlated density is carried on by expanding the correlation function ( p^ in 
terms of Legendre polynomials. This facilitates the calculation in the spherical basis of the s.p. wave 
functions (§). Details are given in reference |l3| . 



4 Applications 
4.1 The input 

We have conducted our study in three different mass regions located around the doubly magic nuclei ^^O, 
""^"Ca and ^'^^Pb. In our calculations we used the same set of s.p. wavefunctions for evaluating both LRC 
and SRC effects. These wavefunctions have been generated by diagonalizing a Woods-Saxon potential 
in a harmonic oscillator basis. This procedure automatically discretizes the continuum. In the RPA 
calculations we used configuration spaces considering all the states below the Fermi surface and three 
major shells above it. We found good convergence of our calculations when multipole excitations up to 
J=10 were used. 
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Table 1: Parameters of the Landau-Migdal force, Vlau and of the zero-range part of the Jiihch-Stony 
Brook interaction Vjs- We used the same normahzation constant Co for both interaction. For ^^O and 
208pb calculations we used Co=386 MeV fm^, while for ^^Ca. Co=212 MeV fm^. The density dependent 
parameters have been changed for each nucleus considered. They are q;=0.50, 0.53, 0.60 fm^^ and Rp— 
3.0, 4.1, 7.22 fm respectively for ^^Q, ^OCa and ^ospb. 

The parameters of the potential we have adopted are taken form the literature ||l^, |l8| and have also 
been used in the FHNC calculations of reference [|l^ , where j-j coupling and isospin dependence of the 
s.p. wave functions were considered. We obtained the charge distribution by folding the pointlike proton 
distributions with the nucleon electromagnetic form factors of reference 

The RPA calculations used to evaluate the LRC effects require an additional input: the particle- hole 
interaction. We have tested the sensitivity of our results to the choice of this input by using two different 
types of interaction. A first one is a Landau-Migdal (LM) contact interaction of the form: 

VLA/(ri,r2) = Co [h{ri) + fa{ri)Ti ■ T2 + go{p)(Ti ■ (T2 + g'o{p)(Ti ■ ct2Ti ■ T2] S^{ri - (18) 

The experience of RPA calculations in the Pb region |]lOj has shown the need of imposing a density 
dependence to the scalar and isospin parameters /o and /g parameters. As suggested in reference [|l^ for 
these two channels we use the expression: 

1 + exp[a{r — Rf)] 

where a and Rp are two free parameters. The factor Co represents the inverse of the density of states at 
the Fermi surface and we consider it as a free parameter used to set the full strength of the force. 

The second RPA interaction is the finite range Jiilich-Stony Brook (JSB) force obtained by 
adding to the zero range piece of Landau-Migdal type, the contribution of one-pion plus one-rho meson 
exchange potentials. This is a finite range interaction and in the discrete excitation spectrum is necessary 
to obtain a good description of the magnetic excitations [^^ . 

The values of the parameters used in our calculations are given in table The parameters of the 
delta interaction are those of reference [|l^ where have been fixed to reproduce the excitation energy of 
the low lying 2+ and 3~ states in ^'^^Pb as well as the position of the the monopole giant resonance. We 
have adjusted the free parameters of the JSB force to reproduce the 3^ state of ^osp]^ 2. 61 MeV. We 
obtained for the B{E3; 0^3) values: 7.536 x 10^ e^fm^ with LM interaction and 7.489 x 10^ e^fm^ with 
JSB, in agreement with the experiment [p3| . 

Even though the parameters of the two interactions have been tuned in the lead region their perfor- 
mances in reproducing the ^^O spectrum are quite good. The low lying 3^ state of ^^O was obtained at 
the energy of 6.95 MeV and at 7.35 MeV respectively for the LM and the JSB interactions. These values 
should be compared with the experimental values of 6.13 MeV. In this case the B{E3;0 3) values 
are 0.97 x 10^ e^fm^ for LM and 1.17 x 10'^ e^fm^ for JSB to be compared with the experimental value, 
1.5 X 10'^ e^fm^ This good result is obtained because the empirical single particle energies are rather 
well reproduced by the spherical Woods-Saxon potential. 

The situation is rather different for ""^Ca. In this case the splitting between the single particle levels 
around the Fermi surface is narrower than the empirical one. For this reason the interaction fixed in 
208pj^ is too attractive. We got an imaginary eigenvalue for the first 3~ state. We reduced the value of 
Co from 386 MeV fm'^ down to 212 MeV fm'^ and we obtained for the 3~ the energies of 1.48 MeV and 
3.20 MeV for LM and the JSB forces respectively. The experimental value is 3.73 MeV. The LM force 
gives a too low B{E3;0 3) value, 1.1762 x 10^ e^fm^ (the experimental one is 1.70 x 10'' e^fm*^ ||2|), 
while the JSB gives 1.22 x 10"* e^fm*^. 
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Figure 2: Differences between IPM charge densities and those obtained considering LRC correlations. 
The full lines have been obtained with the JiiUch-Stony Brook interaction, the dashed lines with the 
Landau-Midgal force. 
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Figure 3: Differences between IPM densities and those obtained considering SRC correlations. The full 
lines have been obtained with the state dependent correlations of reference Q and the long dashed lines 
by using only the central term of this correlation. The gaussian correlation of reference |l9| produces the 
short dashed lines. 
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In figure g we show the differences between the IPM charge distributions and those calculated con- 
sidering the LRC evaluated with the two interactions. We can see that the densities calculated with the 
Landau-Migdal and the Jiilich-Stony Brook interactions are very similar. The differences between the 
results obtained with the two forces are within the accuracy of the data. 

In our model, the input required to evaluate the SCR effects are the correlation functions of eq.(p^. 
In variational calculations these functions, together with the s.p. wave functions, are determined by the 
minimization procedure. In the calculations of reference where the semi-realistic nucleon-nucleon S3 
interaction has been used simple scalar correlation functions were required. In our calculations we 
have considered the gaussian correlations of reference Their expression is: 

/(i)(r) = 1 - Acxp(-Br2) (20) 

and the parameters A and B have been fixed by the minimization procedure. It has been shown in 
reference |l9| ] that, for a fixed nucleon-nucleon interaction, the correlation functions are almost the same 
for every nucleus considered. The obtained values of the two parameters were A=0.7 and B—2.2 fm^^. 
The differences between IPM charge distributions and those calculated with this correlation function are 
given by the short dashed-lines of figure |[ 

We studied the effects of the state dependent terms of the correlation function by using the correlations 
defined by the FHNC calculations of reference In this reference the realistic V8' nucleon-nucleon 
interaction plus the Urbana IX three-nucleon interaction |2) have been used. The difference between IPM 
charge distributions and those obtained with this correlation functions are presented in figure ^ by the 
full lines. We performed calculations switching off the state dependent part of the correlations and we 
used only the scalar part. The results of these calculations are given by the long dashed lines of figure 
^. When the state dependent part of the correlation is switched off we essentially reproduce the curves 
obtained with the gaussian correlation. The contribution of the state dependent correlations diminish 
the difference with the IPM results. These results confirm the findings of references p^ . 

In the calculations we shall discuss in the next sub-sections we shall refer to the LRC calculations 
as those done with the Jiilich-Stony Brook interaction and to the SRC as those done with the full state 
dependent correlation of reference B . 



4.2 Doubly magic nuclei 

The IPM charge distributions of ^^O, '^"Ca and ^°^Pb are compared in figure ^ with those obtained after 
the inclusion of the correlations. The effects produced by both kind of correlations are analogous. The 
charge distribution is lowered in the center of the nucleus and this implies an increase of the rms charge 
radii as it is shown in table 0. 







4^^Ca 




IPM 


2.64 


3.25 


5.46 


SRC 


2.70 


3.31 


5.52 


LRC 


2.98 


3.57 


5.59 



Table 2: Root mean square charge radii in fm. 

The correlations change the occupation probability of the s.p. states with respect to those of the 
IPM. For the LRC we calculated the occupation probabilities by evaluating the expectation values of the 
operator defined in Eq. . The results are shown in figure ^ where we observe that the main changes 
are around the Fermi surface. 

Our LRC correlations produce slightly stronger effects than those found by Gogny |2^ , whose results 
are shown by the white dots in the ^°^Pb panel of figure |. Our LRC effects are strong also in comparison 
with those of the RPA calculations of reference |^ done in the ^"Ca nucleus. These differences are due 
to the different interactions and s.p. bases used. 
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Figure 4: Charge density distributions for the three closed shell nuclei investigated. IPM (full lines), 
LRC (long dashed lines), SRC (short dashed lines). 
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Figure 5: Proton occupation numbers as a function of the single particle energy. The full lines represent 
a Fermi gas distribution. The white dots in the lowest panel show the results of reference pa . 
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Figure 6: Occupation numbers of the valence hole and particle states as a function of the excitation 
multipole. On the abscissa the multipoles are odered for a fixed value of J first with positive and then 
negative parity. The left vertical scales and the upper lines refer to the hole states and the right scales 
associated with the lower lines to the particle states. The dashed lines show the results obtained by 
leaving out the lowest collective 3~ states. 
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In the LRC calculations the change of the occupation probability is due to the coupling of the s.p. 
level with the collective excitation of the nucleus. This is clearly shown in figure ^ where the occupation 
probabilities of the two s.p. levels just below and above the Fermi surface are shown as a function of the 
multipole excitation. It is evident that the occupation probability is strongly modified in correspondence 
to the 3~ states. The dashed lines show the result of a calculation when the lowest collective 3~ state is 
not considered. 

The SRC correlations produce analogous effects on the occupation probabilities of the s.p. levels. 
An evaluation of these modifications is not straightforward as in the case of LRC since in the CBF 
calculations the quantities used are the one- and two-body densities and not the s.p. wave functions. We 
have calculated the one-body densities for the three nuclei considered within our first order correlation 
model by using the multipole decomposition technique described in references jl^, We projected the 
/ multipole term pi(ri,r2) of the one-body density matrix on the s.p. basis by calculating 

nsRC = J dnrl J dr2rlRnijiri)Rnij{r2)piiri,r2) (21) 

The occupation numbers for the s waves obtained with this procedure are compared in table ^ with the 
LRC ones. In our calculation the LRC are more effective in lowering the occupation probability of the 
valence states than the SRC. 
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lsl/2 


0.93 0.95 


0.93 0.94 


0.97 0.93 


2sl/2 




0.78 0.94 


0.92 0.91 


3sl/2 






0.80 0.92 



Table 3: Occupation probability of the s states. 

We obtained for the occupation probability of the 3sl/2 state of ^°*Pb the value of 0.74, which is 
slightly larger than the 0.6 extracted from the (e,e'p) data of ref. even though recent DWBA re- 
analysis obtained a value of 0.7 The nuclear matter studies of SRC and LRC of the occupation 
probability predict values around 0.8 |2|, |2^. As expected our value is lower than this since we also 
include surface vibrations. 

4.3 Semi-magic nuclei 

We have calculated the differences between the charge distributions of isotones differing by one proton 
only, since we expect this quantity to be less dependent from the choice of the s.p. wave functions than 
the full charge distribution. The results of these calculations for the three pairs of isotones considered 
are shown in figure ^ The IPM results correspond to the square of the s.p. wave function of the least 
bound proton state. They are the lpi/2, 1^3/2 and 3si/2 s.p. states for the ^^O, ''"Ca and ^"^Pb nuclei 
respectively. The shape of the curves shown in the figure reproduce the known behaviour of these s.p. 
wave functions. Only the s wave is peaked in the center of the nucleus, the other ones have their maxima 
at larger distances. 

With respect to the IPM results, the correlations diminish the difference in the central part of the 
nucleus and tends to shift the charge in the external regions. In the i^O-^^N and ^°Ca~^^K systems this 
effect produces a shift of the maxima. In the 208p|-,__207rpj gygtcm the inclusion of the correlations lowers 
the difference in the center of the nucleus. As in the case of the doubly magic nuclei short and long range 
correlations produce effects of the same order of magnitude. 

The interest in the calculation of these charge differences is related to the existence of experimental 
data in the lead region. In a single experiment Q the elastic electron scattering cross sections of 
^"^Pb and ^os^pj have been measured. From the analysis of these data the empirical difference between 
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Figure 7: Charge distribution differences between isotones. IPM (full lines), SRC (short dashed lines) 
and LRC (long dashed lines). 
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Figure 8: Differences between ^oep]^ g^^^ 205rpj ^.j^g^j-gg densities. Full lines IPM, dotted lines SRC, dashed 
lines LRC, dashed dotted lines full calculation. In the lower panel the charge densities are multiplied by 
r^ to emphasize the large distance behaviour. The experimental data are from pi pO|. 
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the charge distributions of these two nuclei has been extracted. In the center of the nucleus the IPM 
predicts a larger charge difference than the one deduced by the experimental data. 

A direct comparison between the differences shown in the lowest panel of figure |^ and the experimental 
one is not correct since the ^°^Pb and ^os^pj nuclei have a rich open shell structure which is not considered 
in the ^°*Pb and ^OT^pj system. To describe the two open shell nuclei we adopted the model already used 
in reference The ground state of ^os^j jg described as: 

|1/2+,205t1> = ai(|3si/2 ®|0+,2"^Pb> 
+ a2(|2d3/2 ®|2+,206pb> 

+ a3(|2d5/2 ®|2+,2"6pb> (22) 

where we considered ai=0.86, a2=-0.47, a3=0.30, obtained in reference |^ using first order perturbation 
theory with a 5 interaction. The 0+ and 2+ states of the ^"^Pb are obtained as linear combination of two- 
hole wave functions in the closed neutron core of ^°*Pb |Q. We use the hole- hole amplitudes calculated 
in reference | |3^ within a Tamm-Dankoff Approximation. The ground state amplitudes are rather similar 
to those given in reference which reproduce the empirical values obtained by an analysis of (d,p) 
reactions on ^°^Pb target. The IPM charge difference obtained with these wavefunctions is represented 
by the full lines of figure |^. 

While the LRC have been consistently treated within the open shell structure of the two nuclei p^ , 
this has not be done for the SRC. In this case we supposed their effect to be the same as that calculated 
in closed shell nuclei: 

and analogously for the 205-207 system. With this hypothesis the total charge difference between ^°^Pb 
and 205T1 is: 

where we have indicated with A the difference between the charge densities of the two nuclei indicated 
in the upper index. The total difference evaluated in this way is compared in figure || with the empirical 
one. The inclusion of both long and short range correlations greatly improves the agreement with the 
experiment. In the lower panel we show the differences multiplied by to emphasize the comparison 
at large values of the nuclear radius. Also in this region the inclusion of both correlations improves the 
agreement. 



5 Summary and Conclusions 

We have calculated the effects of long and short-range correlations on the charge distribution of some 
doubly magic nuclei. The SRC have been treated in the framework of the Correlated Basis Function 
theory. We use a first order model in the correlation function. The correlations have been taken from 
FHNC calculations done with realistic interaction. The LRC correlations have been treated within a RPA 
framework. In this case we used an effective interaction whose parameters have been fixed to reproduce 
the ^°*Pb low-energy excitation spectrum. 

We have seen that long- and short- range correlations produce similar effects on the charge distribu- 
tions. In both cases the IPM distributions are lowered in the center of the nucleus and, as a consequence, 
the rms radii increase. In terms of s.p. occupation numbers this effect is produced because both kinds of 
correlations decrease the occupation of the hole states and increase that of the particle states. We have 
shown that in all the nuclei considered the link between s.p. levels and collective low-lying 3~ state is 
one of the major sources of modification of the occupation number. 

The charge differences with nuclei having one proton less have been calculated in order to compare 
our results with the measured empirical difference in the Pb region. Also on the charge difference the 
two kind of correlations produce similar results. Also in this case the correlations tend to diminish the 
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Figure 9: Results obtained with a new set of s. p. states (see text). Upper panel: ^"*Pb charge distributions 
compared with the empirical one |Q represented by the thick line. The thin full line show the IPM result, 
the dashed line has been obtained by summing long and short range correlations. Lower panel: Differences 
between ^°^Pb and ^05^1 charge densities. The meaning of the lines is analogous to that of figure ||. 
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Figure 10: Above: elastic electron scattering cross section on 208p|-, ((jat^ from Below: ratio beteen 

the elastic cross sections on ^os^pj ^^^^ ^^^Pb targets (data from |^). The full lines represent the IPM 
results and the dashed-dottoed ones the results obtained includibg both LRC and SRC. 
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difference in the central region. After a proper consideration of the open shell structure of the ^"^Pb and 
205rpj jj^f-igj^ ^}jg total difference which considers both correlations can reproduce the empirical one. 

Our calculations do not treat LRC and SRC on the same ground. A mixing of ingredients taken 
from fundamental theory, the correlation functions, and effective theory, the RPA correlations, has been 
used in our calculations. In all our work the choice of s.p. wave functions is essential. The sets of s.p. 
states we have used have been taken from the literature and the parameters of the mean field used to 
generate them have been fixed to reproduce at best the s.p. energies of the valence states around the 
Fermi level, and the charge distributions. For this reason we did not make a comparison between our 
charge distributions and the experimental ones. Since the IPM densities are already reproducing rather 
well the empirical ones, the inclusion of correlations worsen the agreement with the experiment. 

It is sufficient a small change in the mean field parameters to obtain a set of s.p. states which can 
reproduce the empirical densities after the inclusion of the correlations. We have done that for ^ospj^ 
where we have changed the Woods-Saxon depth Vq from 60.4 MeV to 62.0 MeV and the radius from 7.46 
fm to 7.4 fm. The results obtained with this new set of s.p. wave functions are shown in figure |^. In the 
upper panel we compare ^''^Pb charge densities with the empirical one [^^ . In the lower panel we show 
that even with the new set of s.p. states the charge density difference between ^"^Pb and ^os^pj^ ^^y^ 
reproduced. 

In the upper panel of figure |l^ we compare with the data of ref. |^ the electron scattering elastic 
cross sections calculated within the Distorted Wave Born Approximation by using the densities of the 
previous figure. The dashed-dotted hue obtained by including both SRC and LRC reproduces better 
the data than the IPM density (full line). In the lower panel we compare our results with the ratio 
of the elastic cross sections measured on ^°^Pb and '^'^^Tl targets [|]. Also in this case the line better 
reproducing the data is the dashed-dotted one obtained by including all the correlation effects. 

We have presented these results to show that we would be able to describe the empirical densities 
with small changes of the mean field parameters. This was not the goal of our work. We simply wanted 
to investigate the relative role played by long and short range correlations, and we found that they should 
be simultaneously considered in order to obtain a successful description of the empirical densities. This 
is our main message. 

The facts triggering our investigation are related to the difficulties found by recent CBF variational 
calculations [Q in reproducing the experimental binding energies and charge distributions of some doubly 
closed shell nuclei. These microscopic calculations, depending only from the realistic nucleon-nucleon 
interaction, provide a good description of the SRC. Our work adds another piece of evidence that the 
adequate treatment of the LRC is the missing ingredient. Hopefully, perturbative corrections, like those 
developed in in references pi, |3^ for nuclear matter, could provide the required description of the LRC. 



Acknow^legments 

We thank P.F. Bortignon and A.M. Lallena for useful discussions. This work has been partially supported 
by MURST through the Progetto di Ricerca di Interesse Nazionale: Fisica teorica del nucleo atomico e 
dei sistemi a molticorpi. 



18 



References 

[1] Cavedon J M et ai 1982 Phys. Rev. Lett. 49 978 

[2] Pudliner B S, Pandharipande V R, Carlson J, Pieper S C and Wiringa R B 1997 Phys. Rev. C 56 
2261 

[3] Wiringa R B, Pieper S C, Carlson J and Pandharipande V R 2000 Phys. Rev. C 62 0144001 

[4] Fabrocini A, Arias de Saavedra F and Co' G 2000 Phys. Rev. C 61 044302 

[5] Mokhtar S R, Co' G and Lallena A M 2000 Phys. Rev. C 62 067304 

[6] Fabrocini A and Co' G 2001 Phys. Rev. C 63 044319 

[7] Rowc D J 1968 Phys. Rev. 175 1283 

[8] Lenske H. and Wambach J 1990 Phys. Lett. B 249 377 

[9] Migdal A B 1967 Theory of finite Fermi systems and appUcations to atomic nuclei (New York: John 
Wiley) 

[10] Spcth J, Werner E and Wild W 1977 Phys. Rep. 33 127 
[11] Co' G and Spcth J 1986 Phys. Rev. Lett. 57 547 
[12] Co' G and Spcth J 1987 Zeit. Phys. A 326 361 
[13] Co' G Nuov. Cim. 1995 A 108 623 

[14] Arias de Saavedra F, Co' G and Renis M M 1997 Phys. Rev. C 55 673 

[15] Clark J W 1979 Prog. Part, and Nucl. Phys. 3 89 

Pandharipande V R and Wiringa R B 1979 Rev. Mod. Phys. 51 821 

S. Rosati 1982 From nuclei to particles, Proc. Int. School E. Fermi, course LXXIX ed. A. Molinari 
(Amsterdam: North HoUand) pp 73-112 

[16] Amaro J E, Lallena A M, Co' G and Fabrocini A 1998 Phys. Rev. C 57 3473 

[17] Rinker G A and Speth J 1978 Nucl. Phys. A 306 360 

[18] Co' G, Lallena A M and Donnelly T W 1987 Nucl. Phys. A 469 684 

[19] Arias dc Saavedra F, Co' G, Fabrocini A and Fantoni S 1996 Nucl. Phys. A 605 359 
[20] Bcrtozzi W, Friar J, Hciscnbcrg J and Ncgclc J W 1972 Phys. Lett. B 41 408 
[21] Spcth J, Klcmt V, Wambach J and Brown G E 1980 Nucl. Phys. A 343 382 
[22] Co' G and Lallena A M 1990 Nucl. Phys. A 510 139 

[23] A. Bohr and B. Mottelson 1975 Nuclear Structure, Vol. II (London: Benjamin) pp 561 and 642 
[24] Afnan I R and Tang Y C 1968 Phys. Rev. 175 1337 

[25] Gogny D 1979 Nuclear Physics with Electromagnetic Interactions, ed H. Arenhovel and D. Drechsel, 
Lecture Notes in Physics, Vol. 108 (Berlin: Springer) p 88 

[26] Quint E N M et aJ 1986 Phys. Rev. Lett. 57 186 
Quint E N M et ai 1987 Phys. Rev. Lett. 58 1088 



19 



[27] Udi'as J M, Sarriguren P, Moya de Guerra E, Garrido E and Caballero J A 1993 Phys. Rev. C 48, 
2731 

[28] Pandharipande V R, Papanicolas C N and Wmbach J 1984 Phys. Rev. Lett. 53, 1133 

[29] Benhar O, Fabrocini A and Fantoni S 1990 Phys. Rev. C 41, R24 

[30] Frois B et al 1983 Nucl. Phys. A A396 409c 

[31] Zamick L, Klemt V and Speth J 1985 Nucl. Phys. A 245 365 

[32] Klemt V and Speth J 1976 Zeit. Phys. A 278 59 

[33] Kuo T T S and Herling G H 1971 JVavaJ Research Laboratory Memorandum Report 2258 

[34] De Jager C W and De Vries C 1987 At. Data Nucl. Data Tables 36 495 

[35] Heisenberg J et al 1969 Phys. Rev. Lett. 23 1402 
Eutener M et al 1976 Phys. Rev. Lett. 36 129 
Frois B et al 1977 Phys. Rev. Lett. 38 152 

[36] Fantoni S and Pandharipande V R Nucl. Phys. A 427 473 

[37] Benhar O, Fabrocini A and Fantoni S 1989 Nucl. Phys. A 505 267 



20 



